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ABSTRACT 


This thesis addresses the problem of determining the 
optimal number of spares for a multi-item inventory system 
with a procurement budget constraint. Various inventory 
models are considered with objective functions like time- 
weighted units short, units short, essentiality-weighted 
units short and pseudo-availability. Solution algorithms 
are derived using the generalized Lagrange multiplier 
approach and a marginal analysis approach. 

Sample data and output results are provided and compari- 
sons of the alternative models are given. Finally, a dis- 
cussion and example is given:of the use of the models as 
a means of estimating the budget required to attain a speci- 


fied level of performance. 
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Ef. IJINPRODUCTION 


In today's world, while all systems are becoming more 
and more sophisticated, the control and maintenance of 
inventories of these systems 1S a problem common to all 
enterprises and military services. In private and commer- 
cial concerns the effective control of inventories can re- 
sult in decreased costs, increased sales and profits and 
consumer satisfaction. In the military proper management of 
inventories may contribute to increased availability and 
readiness, decreased inventory investment and system costs. 

For each component of each weapon system two fundamental 
guestions must be answered: 

(1) When to replenish the inventory; 

(2) How much to buy for the replenishment. 
In order to answer these questions, many inventory models 
have been developed in the past 30 years. See for example, 
Hadley and Whitin [Ref. 1], Muckstadt [Ref. 2] and Eriksson 
(Ref. 3]. Most previous work solves a variety of cost 
minimization problems considering expected values of steady 
state variable costs associated with shortage cost, ordering 
cost and storage cost. 

Such models may be appropriate for the commercial sector, 


but are not always appropriate in the military world. 





In the commercial sector, the objective function of the 
inventory model is to maximize profit or minimize the average 
annual costs. Non-cost oriented objective functions frequently 
are used in the military inventory systems. For example, 
attempts are often made to maximize availability or fill 
rate, Or minimize the number of backorders or expected time 
weighted stockouts, or minimize the probability of a stock- 
out with a budget constraint. 

Obviously costs are important in every inventory model. 
However, Many real-world inventory problems are so compli- 
cated, one cannot represent accurately the real situation. 
Thus, some simplifications and approximations are used when 
constructing a mathematical model of any real world system. 

If this is not done, the results obtained by use of the model 
can easily lead to operating rules which are worse than those 
currently in use, worse than those which could be derived 
from simple heuristic intuitive considerations. 

Many of the inventory problems are viewed as single period 
problems. For example, initial provisioning, allowance list 
determination and the fly-away kit problem are single period 
problems. These models are perhaps the simplest of the models 
in which demand is treated as a stochastic variable. 

Reasonable objective functions in these models are to 
maximize performance subject to a constraint on the resources. 
Typical measures of performance might be availability, time- 
weighted units short, fill rate, the number of backorders, 


and mean supply response time. 





This thesis considers various single period models which 
attempt to maximize performance subject to budget constraint. 
Chapter II describes the general single period problem 
and introduces the method used in this thesis of solving 

those problems. 

Chapters III and IV develop the time-weighted units short 
model and the availability model, and explain the solution 
procedure. Sample data runs for both models are provided. 

Chapter V provides a comparison of models considered in 
the thesis and discusses some of the properties of each 
model, and Chapter VI discusses the use of models for purposes 
of determining the budget required. 

Chapter VII summarizes the results of the research and 


concludes with some suggestions for additional research. 
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II. THE GENERAL PROBLEM 


In this chapter, we consider the general single period 
modei as a process for transforming resources into new dis- 
tributions of inventory positions over the line items in 
the inventory. 

The essential problems of control in a line item-inventory 
control system with multiple line items are: 

(1) How much resources to commit at a point in time; 
(2) How shall these resources be allocated to achieve 
system objectives. 

In a typical continuous review inventory system, we can 
determine the optimal order quantity (Q) and reorder point 
(x) for a given item by minimizing the average annual variable 
costs. But in applying this theory to the real world inven- 
tory systems which consist of multiple line items, it is 
frequently the case that resulting minimum cost solutions are 
not feasible because of a budget limitation or some other 
constraint. Thus, in a constrained multi-item inventory 
system, the typical continuous-review policy is sometimes 
inappropriate. In the following section we discuss several 
objective functions to guide the line item inventory control 
system in determining how to allocate availabile procurement 


funds at a particular replenishment epoch. 


le): 





A. GENERAL FORM OF OBJECTIVES WITH CONSTRAINTS 

Consider the case in which an administrator, responsible 
for the replenishment decisions, determines replenishment of 
stocks of various line items on a periodic basis. Suppose 
that a fixed amount of procurement budget has been allocated 
to the replenishment epoch at hand and that a target number 
of reorder actions has been established as a working con- 
straint for the allocation epoch. The administrator's task 
1s to transform the available resources into replenishment 
orders for different items. 

1. Measure of Effectiveness and Objective Function 

Daeschner [Ref. 5] examined the constrained line-item 

allocation problems. He considered several possible objective 
functions which can be adapted to the case where unsatisfied 
demands are backordered and to the case where unsatisfied 
demands are lost sales. Let [ > 0 be the penalty (reward) 
per unit for item j and let ve be the demand for item j ina 
period. Let x. be the inventory position for item j after 
ordering ina period. Let oe = rT Then the number of sales 


for item j in the period is given by 


The expected sales for item j is therefore 


12 





x. 
Gyaehee—- ad.) + X. ptD. > X. 
ares ‘ln Zia ag 
j 


which is equivalent to 


We assume that the inventory system seeks to minimize the 
expected penalty incurred, or, equivalently, to maximize the 
expected penalty avoided. Mathematically, the objective is 


to maximize 


N 00 
Z(X) = } w.(E(D.) - (dl. seek, pig sds) 
7 j=1 J J ache — sa) 


Several interpretations and uses of the penalty coefficient 
ue are possible. Four are illustrated in Table I below. Each 
reflects a formulation of system objectives which has been 
adopted or considered by the Navy Supply System. Daeschner 
(Ref. 5] also considered ue as a linear combination of various 
coefficients in his line item allocation model. 

There are many other types of objective functions which 
are currently used in the military. For example, 


(1) Minimize units Short in a given period 


N oO 
Z(X) = eee pi oe =.) 
a Py q=K,+1 } J J J 


iis 





TABLE I 


INTERPRETATIONS AND USES OF T 


renakty Coemracients Objective 
ue = fs Maximize expected sales from 
stock. 
ae = i. Maximize the expected requi- 
Sitions filled Whe = average 


quantity of item j demanded 


per requisition). 


Tn. = 1 Maximize the expected number 


of units issued from stock. 


T. = ae + TMNIS Maximize expected customer 

- TMISS waiting time per unit avoided 
by issue from stock, where 
ae is the lead time for item 
j, TMNIS is the calendar time 
anticipated to process a 
request and TMISS is the time 
to affect issue from stock of 


a demand, available item. 


(2) Minimize time-weighted units short 


Z(X) = 5 TWUS, (X;) - E, 
j=1 
where; 
TWUS (X, ) = time weighted units short 
E. = essentiality 


aL 
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(3) Maximize system availability 


N 
Z(X) = I A. (X,) 
i=l 
where: 
A. = Availability for each item. 


i 


The objective functions (2) and (3) will be explained in 
Chapters III and IV. 

2. The Line-Item Allocation Model and Solution Procedure 

In the previous sections, many kinds of objective 

functions are introduced. If we define them correctly, those 
objective functions can be solved by various techniques. It 
is evident that an actual inventory system with limited re- 
sources might be unable to carry out a prescribed inventory 
policy if either the amount of procurement funds available 
or the number of replenishment actions exceed the available 
resources. The problem is made more complicated by the fact 
that the objective functions are "non-linear" and the require- 
ment that the a must be integers. The problem is stated 


mathematically as 


max Z(s) 
N 
(Al) Se). tz chs & 
j=1 J 
N 
) H(s.) < R 
Sere 


ILS) 





s = (S,1Sor---/N): integers 


=e = the unit price of item j 
ae = the number of buys of item j 
H(s.) = 1 ie gSie> 0 

J 3 


0 otherwise 


B = the procurement budget limit at the realloca- 
tion epoch 


R = the maximum number of individual procurement 
activities allowed in the present allocation. 
To solve the problem (Al), the generalized Lagrange multiplier 
(GLM) method of Everett [Ref. 4] can be used. Using this 
method, the problem can be reexpressed as 


N 
(A2) max L(S,A) = 2Z(s) - A, (( ) c.S.) -B) 
s — = ee 


N 
= AQ (CD 


H(s.)) - R) 
j 3 


a 


A. > 0 with optimal solution S*(A). 


Se sane Aye 5 


Problem (A2) is the Lagrangian problem associated with (Al). 
Using Everett's theorem, one can determine a bound on the 


optimal solution, Z(s*) to be 


Z(s*) < Z2(s*(A)) = A, (B(A) -B) - AZ(R(A) -R) 
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where 


N 

Beoue= jf Cese(X) 
a 
N 

Ri) =) HS. (A)) 
| 


ances *(\) sisethe optimal,solution vector for a given pair 


(rae, Xr In solving problem (A2), we can separate the N- 


1’ 9) ° 
variable optimization problem into N one-variable problems. 


Choosing trial values of hy and ho we maximize 


A3 Tix ciel = 2Z. -) - A.C.s. - AWAH(S. 
(A3) oS r) 7 tsq) = 2 (s5) 


When considering the integer nature of decision variables, 
the optimal solutions for (A3) are determined by finding the 


values s such that 
Meee y=] te. ,A) < 0 and L.(s.,A) - L.(s"-l,A) > O. 
, (S#+1,.) ,(S3-N) < 5 (85-0) ,(S4-1/A) 
Thus S; is the smallest value such that 
AL.(S.,A = L.(s.¢+1,A) - L.(s.,dA <8) 
5 (857A) ,(83+1,A) 5(S;/A) 


In order to get an optimal solution, Daeshner [Ref. 5] used 


an interactive computer program, which evaluates the current 
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optimal solutions with hy and Ao - Each time, the user can 


select a pair A ho and objective function type to be con- 


1? 
Sidered. Then the user is provided with output which indi- 
cates the budget consumed, the number of stock replenishments 
generated, the achieved objective function value and a maximum 
attainable value for the objective function. 

After examining the output, the user can modify the 
input parameters and continue or terminate the run. Decreas- 
ing the non-negative multiplier values tends to use more of 
the corresponding resources, increasing the values used, less. 
When the replenishment actions generated by a pair of values 
(A, A>) exactly consume the available resources, B and R, 
the solution is optimal. Frequently exact equality may be 
impossible because of integer nature of the problem. Thus 


the solution obtained may not be optimal, but the difference 


is not likely to be significant. 


B. AUTOMATING SEARCH ON THE LAGRANGE MULTIPLIER 

The interactive search method cannot guarantee an optimal 
solution and it requires trial and error to get the approxi- 
mate optimal solution. Consider the case in which there is 
only a single constraint, with the same type of objective 


functions. The mathematical program is then 


max Z(s) 


(Bl) 
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no) 


(S)7So7-++ Sy) = integer number of buys 


[9 
I 


budget limit 


Cc. = price of item i. 


We can rewrite the above equation using a Lagrange multiplier, 


ass 


(1) Iiey0). = Zs eaner 
7 


[|r 


c.s. -— B] 
1:21 


Then separate the equation. 


N 
= a (Z(s.) - Oc.s.) + 6B 


(2) L(s.,,s 
1 al 


grees rSy) 


Equation (2) can be maximized by maximizing each sub- 
objective function. If Z(s,) is differentiable with respect 


to each Si) the optimal solution is obtained by 








adZ(s.) 
dL _ tet a 
S.c5) ds. Oc., a 4 | 
1 ale 
Thus set x = 0 and get 
a 
dz, (s.) 
: a ds. /C; 
a 
N 
a 


that @ can also be interpreted as a shadow price for the 


objective function: i.e., @ = 3Z2/3B. Due to the integer 


ih 





nature of Si it is often impossible to get an exact optimal 
solution. Difference equations must be used because the 
region in which the solution is desired consists of a set 


of discrete points. Therefore, let 


3) AL, (s; ,8) = L.(s;,8) = L,(s,-1,6) 
= = = = era 
Z.(s;) OC .S5 Z.(s, 1) + 8c. (s, ts) 
= AZ, (s,) = c, 8 


We know that Equation (3) is a concave function at the point 
Eo. The optimal selution must satisfy AL, (s;,8) > 0 and 
AL, (s;+1,8) < 0. Thus the optimal solutions are given by 


finding the largest s,'s such that 
AL, (s;,6) Eso 
or equivalently 
(4) AZ. (s;) — c, 6 > 0 eile, Sc: Nc 


The Lagrangian multiplier @ can be found by the following 


search algorithm. 
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STEP 1. Find an initial upper bound P ae Let all S, be 
assigned zero at the beginning and find the 
change of objective function per unit dollar as 


a result of increasing to one unit. 


AZ, (1) 
1 Cy 
AZ, (1) 
8. 7 Cc 
2 
AZ (1) 
a Nn 
— 
N 


where AZ, (1) 


Z. (4) ze Z, (0). 


Because of decreasing marginal returns or objective function 
values and because of the interpretation of 98, an upper bound 
on @ is given by: Oe, = max{6, ,0 


gr-++79,]. 


STEP 2. The initial a will be 


where oF = 0. 
Find for each i, the largest S; so that 


AZ, (s;) 


Cc. 
a 


| 
O 


and evaluate the objective function and the budget 
required. 
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STEP 3. If the budget used is greater than the given 


budget, let 


(@ + 8 
o — 2 
otherwise 
— eo + 8, 
i 2 


Each time update the S vector, the objective 
function values, the upper bound of the objective, and the 


amount of budget consumed. 


STEP 4. Stopping rule. 
Stop when the used budget is equal to the given 
budget or the difference between the current upper 
bound and the objective function value is less 
than some limit (¢€). Otherwise go to step 2, 
and continue until the above conditions are 
satisfied. 


A FORTRAN program for this algorithm is given in Appendix B. 


C. MARGINAL ANALYSIS PROCESS 

The theory of marginal analysis has been used in many 
inventory models when resource constraints are active. [In 
an economic sense, AZ. (s,)/c; can be interpreted as the 


marginal increase in the objective function per dollar spent 
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achieved by adding one more unit of stock. It is reasonable 
for an inventory controller who has a scarce resource such as 
a procurement budget to buy an item which gives the maximum 
benefit per dollar spent. 

By using a simple computerized algorithm, the line item 
allocation problem can be solved easily. The first step is 


to set all Ss, = O and compute 


= a AZ, (s+) | 
’ peee ge i 
i ma 2 <n 
If the maximum is taken on for item j, set S. = 1 and 


deduct the unit price for unit j from the budget. The second 


step is then to recompute AZ. and then find 


Agets.) A2.:(s.) 
1 1 } 


—— = g, 
The next unit is assigned to the index j where the maximum 
is taken on, etc. This is continued until adding an additional 
unit exceeds the budget constraint. It should be noted, how- 
ever, that the method described does not insure optimality 
(Ref. 1]. Specifically the method may stop too soon. If the 
item i selected from the marginal analysis has a Cc. value 
greater than the remaining budget, the procedure terminates 
even though some other item j may have a oe value less than 
the remaining budget. An obvious improvement in this area 


could be the inclusion of a subroutine that would select 
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from the remaining items the best one from those having oe 
smaller than the remaining budget. A FORTRAN program for 


performing this marginal analysis is provided in Appendix C. 


D. SAMPLE DATA RUNS OF UNITS SHORT MODEL 

A weapon system consists of 10 components. The system 
manager wants ;to minimize the number of units short by supply- 
ing spare parts to Support the weapon system. Suppose that 
the demand rate, lead time, price and essentiality code for 


each item 1 are known. The objective function can be expressed 


by 
er 
Minimize Z(s) = (d--s.)p(D.=d.)E. 
feu) z i=l d.=s.+1 7 ? a) +e 
fe ese 
10 
subject to 2 s.c, < B 
i=l 
where: 


E; = essentiality code 
B = budget limit 
“A;T; ae 
e er. te) - 
p(D.=d. ) = ee 
ele (d.)! 


The approximate solution of (Cl) can be obtained by the marginal 
analysis method. Table II shows the computational results 


for this system with known input data. 


24 





TABLE If 


THE RESULT OF UNITS SHORT MODEL 


No. Ay Lead Price Essen. Alloca- Units 
Time (S$) tion Short 

il Dhow 0 16 146%, 2 0 340 0.0007 

2 Oreck gO 208.0 dO 140 0.0048 

3 SO iO 100.0 1230 2.0 1.2489 

4 25. 0 150 2.0 370 42.0 0 0013 

5 0 Leo a0 0 a0 0.0007 

6 OS 1.0 a0 3730 4.0 0.0002 

7 SLC) ace) Ae 0 10 1.0 20 OF 00L2 

8 5.0 lee 100.0 1.0 4.0 1.4368 

9 120 6 50.0 1.0 370 O02.33 

10 2.0 ZO 100.0 0 Za 0.5413 


Table II shows several properties of the units short 
model. First of all, more than one unit short in a year 
occurs in the high cost items (items 3 and 8). Second, low 
demands and low price items are allocated enough. Items 2, 
5 and 6 are allocated more than five times their mean demand. 
Also this model tends to stock more of the high demands and 
low price items. 

Finally, the essentiality weights cause greater alloca- 
tions to be provided to those items with high essentiality 


than would be provided with equal weights. 
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Other results include: 


Total objective value 520 
Shadow price 0.001899 
Budget limit $1170 
Budget left SOO 
NICS. Ee. 
The shadow price is the last maximum value of = 
a 


It can be interpreted approximately as the amount of 
decrease in the objective function achieved by adding one 


more dollar. 
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IIIT. TIME WEIGHTED UNITS SHORT MODEL 


A. DESCRIPTION OF MODEL 
In the previous chapter we have discussed various objec- 
tive functions and solution methods for the single period 
inventory problem. In the military services, many measures 
of effectiveness have been used to indicate system performance. 
Among these measures are fill rate, availability, mean supply 
response time, the number of stockouts, and time-weighted 
units short. In this chapter we consider a model which 
minimizes time-weighted-units-short (TWUS). 
Suppose that a weapon system consists of n components 
and the objective is to allocate a given budget for spare 
parts so as to minimize time-weighted-units short for the 
entire system. Assume that 
(1) procurement lead time and repair lead time are known 
constants. 
(2) demands for each installed unit have a known 
distribution. 
(3) the total amount of procurement budget available to 
spend on all components is fixed. 
(4) the objective is to minimize essentiality weighted 
TWUS. Mathematically, the model can be written 


as 


ou 





n 
min ) fTWUS.(s.)E.-: 
: ae 


_ SLT 
n 
Set . » cs, < B 
i=l 
where: 
TWUS , (s;) = time weighted units short when there 
are s, units for item 1 
SLT = total sum of lead time demand 
n 
(fer) 
ol) ss 
E. = essentiality code for item i 
B = budget limit ina given period 
Cc; = price of each item. 


In the above problem, if the TWUS is properly defined, this 
model will be solved easily by using the methods explained 


in Chapter II. 


B. POISSON DEMAND CASE 

We shall now determine an exact expression for the 
TWUS ; (s,) for the case in which demands are Poisson distributed. 
Let the mean rate of demand be hs units per year and the 
lead time be a constant T;. In addition to treating the 
demand variable as being discrete, the number of buys S; 
also will be treated as a discrete variable. Thus if D; 1s 


the lead time demand item j: 
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e te i 
A (s.)1 = 
L 
i p(s,;A,T;) 
Let 
p(s,) = prob(D,>s,) = ; L p(d,;A.T,) (2) 
i "a 


If there are S5 units of stock for item i, Richards and 
McMasters [Ref. 8] show that the expected time-weighted 


units short in (0,T,) is given by 


S,(s,+1) 


Be 
ee = ee 
E[TWUS (s,)] = xip(s,+1) [A,T, 2s, + Wh J 


For those cases where the expected lead time demand is 
large, the Poisson probabilities in (3) can be approximated 


by a normal distribution with mean A,T; and variance 


go? = A\.T.. Let 
1 i. 2 


g(x) = exp (-x*/2) 


— 
V2T 
be the standard normal probability density function and let 


d(x) = f $(u)du be the complementary cumulative distribution 
x 
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function for the standard normal. Then expression (3) can 


Memrewritten in terms of the normal probability function as 


follows: 
T, s,+1-\,T, s.(s.+1l) 
E{TWUS,(s,)] = == [AT -2s, + TT ] 
a Pe 
1 S,-A;T; 
ez a ee aa)? (4) 


This expression should be used in those cases in which A;T, 
is large. We have developed the expression for the expected 
time-weighted units short ina period of length T. when there 
are s. units of stock for item i. 

In the next section we write the expression for the total 
essentiality-weighted time weighted units short over all items, 
and we provide a solution procedure for allocating the given 


budget optimally. 


C. SOLUTION PROCEDURE 
The mathematical program for the time-weighted-units short 


problem is: 





l N E,T; _ s,(s,+1) 
(Cl) min Z(s) = sit 4, 5 ON re i a 
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To solve this problem we can use the Lagrangian multiplier 


technique. Let 


N 
7 Z,(s;) ~ pf. c.s,) G5) 


rr] 


L(S)/S57-+++S,78) = 
a 
Here Equation (5) is separable in the items, and minimization 
of the total objective function is accomplished by minimizing 
the individual functions Z;(s;) subject to budget constraints. 


Consider a single item i. Let 


AL; (s;) = L;(s,;-1) = L(s,) 
= Z,(s,-1) -2,(s;) +8c;s,; - 8c; (s;-1) 
= AZ,(s,) + 8c; (6) 
where 
AZ.(s;) = E,(Twus,(s,-1) - Twus,(s,)] (7) 


As shown earlier 


TWUS (s-1) 5{P(s) [AT-2 (s-1) +3452) ] 4p (5-1; 47) (7T$641) } 


2 
eS 7 $s -s ES) 7 - 
= s{P(s) [AT 2s+= +2] + sap(szAT) OAT st+1) } 


aa. 





TWUS(s) = ={P(s+1) [AT-25+S‘S41) + p(S;(T) (AT=s) } 
_ Ts Reese peesiices)) 
= x{P(s) Px 2st—Sa]) p(s;AT) CAT 2st) 
+ p(s;AT) (AT=-s) } 
so that 


TWUS (s-1) - TWUS(s) 


a 2 ak 
= x{P(s) [2 += 3] +p(s; AT) [s - a +yq+AT-2s 


2 
Ss +s 


7" AT + sj} 


+4 





= 2 2 
= 5(P(s) (2 - $5] + $= p(s;AT)] 


= P(s) (T - 5) + > p(s;AT) (8) 
Substitute Equations (7) and (8) into (6). Then 
S. Ss. 


AL, (s,) = BE, ((T, - Tae sh. iT; i) ty, P (sy ;A5T5)] 


+ Oc, (9) 
The optimum solution Ss; is the largest Ss. such that 


AL.(s,) 2 0 
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or equivalently, 


ga 8841) 2; '8;) 7 —— 
i i i 
Ss. 
+ AG p(s,)] > =e 


The basic algorithm for solving this problem was explained 
in the previous chapter. A computer program for searching 


for 8 is provided in Appendix B. 


D. SAMPLE DATA RUNS 

Consider a weapon system which consists of 10 components. 
Suppose that the demand for each component is Poisson dis- 
tributed with parameter h; and the lead time is known constant 
T,. Let the budget available for procurement be $19224. 
Table III shows the optimal allocations provided by the TWUS 
model. 

The allocation given when the demand distribution is 
approximated by the normal distribution is also provided 
in Table III for comparison. (For comparability, the variance 
for the normal distribution is taken to be the same as the 
mean). Comparing the results, we observe that the normal 
case buys more of the high demands low cost items. There 1s 
a small difference in the allocation for items 8 and 9 which 
are more expensive than the others. For demand rates less 


than 10, the usefulness of the approximation is questionable. 
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Teds, LLL 


OPTIMAL ALLOCATION FOR TWUS MODEL 


Item d. T. Cc. (S) E. Poisson Normal 
1 10 1 10 = ey 20 
2 100 1 20 cei) 13 120 
3 1S 1 80 zt 12 cue 
4 20 1 2 i. A 36 
5 50 1 5 uF 65 aa 
6 80 dh 30 die 90 96 
' d 20 1 1 1 So aii 
8 5 1 200 i. 17 Sy 
9 75 2 100 Hl 77 74 
10 10 iz ia 1 14 16 


The resulting values of the objective function for the 


optimal solutions are: 


Poisson Case Normal Case 
eis *) 0.00094 0.0015 
Shadow price (@*) 0.00015 0.00055 
Budget limit $19224 $19224 
Budget left 0 0 


The objective function value for the Poisson demand case is 
less than the normal demand case. The main reason for the 


difference is due to items 8 and 9. 
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IV. PSEUDO AVAILABILITY MODEL 


A. DESCRIPTION OF MODEL 
In the previous chapter, the TWUS objective function was 

introduced as a means for allocating a limited budget. 
Operational availability is a widely stated measure of the 
operational readiness of military forces and weapon systems. 
Thus, it 1s appropriate to consider stockage models with an 
availability objective as a means of allocating limited 
resources. 

The most direct and meaningful measure of the influence 

of peace time operating stocks on readiness is weapon 

system (or end item) availability. We use the terms 

availability, end-item availability, and weapon system 

availability interchangeably to mean the probability 

that an end item, such as a tank or an aircraft, 

selected at random, is not waiting for a component 

to be repaired or shipped to it. [Ref. 6] 
Many authors have attempted to determine stockage levels for 
components by maximizing equipment operational availability, 
subject to a budget constraint. See, for example, Jee [Ref. 


7]. Usually, the availability for component i is defined by 


the ratio 


MTBF. 
A; = MTBF. + MITR. + MSRT.(s.) 
1 wl ao 1 
where: 
MTBF . = the Mean Time Between Failure of item 1; 


a5 





MTTR, the Mean Time To Repair for item i; 


the Mean Supply Response Time for item 1. 


MSRT.(s.) 
lel 


If the weapon system is assumed to consist of the n components 
all arranged in series, then the system availability is the 
product of the individual item availabilities. (This assump- 
tion means that the system will fail if any of the components 
fails.) With this assumption, the allocation problem, stated 


in terms of system availability is: 


n 
(Pl) max fea cs. ) 
n 
Sse. pg) casi <  B eee. 4. s., N 
; ii - 
i=1 
where: 
A, (s,) = the availability of item i having S; 
units in stock; 
B = the budget limit; 
Cc; = the price for each item il, 


In the expression of A.(s.), the term MTBF . is the recipro- 
cal of the failure rate 1, MITR, is assumed to be independent 
of the decision variables and the available funds and MSRT, (s;) 
can be expressed in terms of TWUS (s,) as 

a 


MSRT,(s,) = 3 Pe Se 
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Thus, the main determination of availability from the point 
of view of the supply system is MSRT.(s,). Many techniques 
for solving this model have been developed. In the next 
section we represent an algorithm for solving the availa- 


bility model by using the marginal analysis method. 


B. SOLUTION PROCEDURE 

The model (Pl) 1s not additive in the individual com- 
ponent availabilities but is converted into an additive 
function by transforming the objective function. Taking the 
natural log of the objective function, the model can be ex- 


pressed in the following way. 


n 
(P2) max ) ny Ag ts.. ) 
S i=l ae, 
n 
sate ) c.Ss, < B ieee... N 
ia 


Now the model (P2) is separable for all 1 and maximization 
of (P2) yields the same solution as maximization of (Pl). 
The marginal analysis method selects an item which gives at 
each step the greatest increase in log(A,; (s;)) per dollar 
epent. 


STEP 1. Start with zero units for all items. 


STEP 2. Compute the increase in log availability per dollar 
spent as a result of purchasing one additional 


Uiaec < 


a7 





Seen ait z 
; 7 el” A, (s,) Qn A, (s,)] 


MOBE . 
al 


in Aj (sj) = inliggyFMPTR, + NSRT|(s,)- 


STEP 3. Select that item i corresponding to the maximum 


ratio. 

Atn A.(s,) Afn A.(s.) A&n A (s_) 
eee 
all s, i 2 n 


STEP 4. Increase the number of units stocked for the item 
selected at step 3 by one additional unit if the 
unit price is less than the amount of budget 


remaining. 


STEP 5. Update the S vector, the MSRT(s) expression and 
decrement the available budget. If the remaining 
budget is greater than the cost of the cheapest 


item, Go to Step 3. Otherwise, Stop. 


In the following section, we will illustrate this procedure 
with a sample system. The computer program is provided in 


Appendix C. 
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C. A NUMERIC EXAMPLE FOR THE MSRT MODEL AND THE AVAILABILITY 
MODEL 


In the expression for availability, the MTBF and MTTR 
terms are not functions of the number of spare parts. There- 
fore it is commonly believed that maximization of system 
availability is equivalent to minimization of mean supply 
response time. However, this is not the case, as shown 
below. 

Suppose a weapon system consists of three components and 
the demands are Poisson distributed with parameters Aye do 


and i respectively. The lead time is a known constant and 


er 
the components have essentiality codes E.- The unit price 
and MTTR are known and the budget is limited to 20 dollars. 


This information is summarized in Table IV. 


TABLE IV 


INPUT DATA FOR EXAMPLE 


ITEM ds C5 MTTR, E; T; 
a i: 5 0.0274 ib So 
Z Ce 5 OU 0zy 3 Eas 
3 10 1 0.0054 i 0 


To solve MSRT minimization problems, we first determine 
the MSRT's for all possible cases. These values are provided 


in Table V. 
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TABLE V 


MSRT DATA FOR ALL FEASIBLE SOLUTIONS Ss 


MSRT (s, ) ITEM 1 ITEM 2 ITEM 3 
MSRT (0) 0.9482 0.9837 0.6 

MSRT (1) 0.3161 0.1967 0.5 

MSRT (2) 0.0708 0.02 0.4099 
MSRT (3) 0.0132 0.00227 0.3298 
MSRT (4) 0.0021 0.0002 0.2596 
MS RT (5) 0.0003 0.1992 
MSRT (6) 0.1485 
MSRT (7) 0.1072 
MS RT (8) 0.0746 
MSRT (9) 0.0500 
MSRT (10) O20322 
MSRT (11) 0.0199 


Using the solution procedure described in the previous chap- 
ter we determine the optimal solution to be as shown in Table 


VI. 
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TABLE VI 


THE ALLOCATION OF SPARE PARTS FOR MSRT MODEL 


ITEM 1 2 3 
USED 
ALLOCATION Adz, (s,) AZ, (s5) AZ,(s,) BUDGET (S$) 

co, 0,0) 0.12462 0.47216 O.1 0 
Ae. ,1.,.0) 0.12462 0.1040 Ouse]. 5 
fie, ,0) 0.04905 0.1040 Oia. 10 
1, 2,0) 0.04905 0.01265 Oi1. 15 
fen 2 1) 0.04905 0.01265 0.09 16 
fle, 2,2) 0.04905 0.01265 0.08012 17 
wi, 2, 3) 0.04905 9.01265 9.07022 18 
(1,2,4) 0.04905 0012'S 5S 0.06039 19 
ie 5) 0.04905 0.01265 6.0507 20 

The optimal solution for MSRT model is (1,2,5). Repeating 


the analysis for the availability objective function we 
obtain the results provided in Table VII from the marginal 


analysis procedure. 


ug 


AZ, (s;) 
i 


di 
ae A,(s,+1) - &n A, (s,)] 
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TABLE VII 


THE ALLOCATION OF SPARE PARTS FOR AVAILABILITY MODEL 


ITEM 1 2 3 
USED 
ALLOCATION AZ, (sy) AZ, (So) AZ,(s,) BUDGET 

ee, 0,0) Orr 7 LZ 0.1294 wei 29 0 
cor 0, 1) Om077 12 0.1294 OF Vell it 
(om0,2) Ore 7 LZ 0.1294 0.1689 Z 
wo, 0 , 3) Oe Oy 72 0.1294 Cry 59 3 
(0,0,4) Or. Cerne 0.1294 0.1808 4 
0,5) OrO77 2 0.1294 0. Lge 5 
c0.,.0 ,6) ORO 7712 0.1294 Om 77 6 
0,7) Oro 702 0.1294 0.1660 a 
(0,0,8) O07 712 0.1294 0.1469 8 
(0,0 ,9) 020772 0.1294 Orel 217 9 
wor, 1,9) Omen? 12 0.0447 Oe Ly 14 
mo, 1,10) Dec77i2 0.0447 0.0938 iS 
co, 1,11) Owe 2 0.0447 0.06702 16 
Or, 7,15) 20 


The optimal solution for the availability model is (0,1,15). 
Comparing the results of the two models, we see that the 
availability model allocates more units to the high demand 


lower cost items than the MSRT model. 
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D. SAMPLE DATA RUNS 

Suppose that a weapon system consists of 10 components 
and the demand of each component is Poisson distributed with 
parameter Aye and lead time To, mean time to repair MTTR, 
are known constants. In order to maximize the availability 
of spare parts with budget constraint, we can use the modi-~ 
fied Availability model (P2) instead of (Pl). By using the 
computer program in Appendix C, this problem can be solved. 
Table VIII provides the allocations of spare parts in the 


Availability model when the budget is 1170 dollars. 


TABLE VIII 


THE ALLOCATION OF SPARE PARTS FOR THE 
AVAILABILITY MODEL 





ITEM =A, 6 6T, cy (CS) E, MTTR, ALLOCATION A, (s,) 
io 1.0 1.0 TORO eels 0 860510137 4.0 0.986 
ce i0md, 1.0 B0.0 io 020274 20 0.997 
fmese0) 61200606 108.006=CUCd SC 37 3720 ons2i 
a 25.0 1.0 2208 3.0. 020622 a7 30 0.327 
a 170 1.0 50m, 10 § “Onoa74 50 0.973 
pe 0.5 1.0 5-0 3.0 0.0027 4.0 0.999 
7 1070 1.0 120) 170 050054 PaO 0.949 
ae 520 150 #410050 1:50 #©«%0,0411 sO 0.538 
9 1.0 1.0 50.0 1.0 0.0082 BO 0.987 

uo meamO 1.0 4100.0 1.0 #£«0.1370 20 0.697 
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From the table, one can see that the availability of an 
item is greatly influenced by the MTTR term (see item 4). 
The availability for that item never exceeds 0.333 even if 
the MSRT is zero. We also observe that the availability 
model tends to stock the high demand low cost items. 


The objective function for the optimal solution is given 


oy: 
Total obejctive value 0.08999 
Shadow price 0.000261 
Budget limit $1170 
Budget left $0.0 


A comparison of the above results with the allocation given 
in Table VIII shows that the total availability is relatively 
low even though most of the items have high availabilities. 
Also as mentioned above, when the MTTR data for an item is 
large relative to the MTBF, a high availability cannot be 


achieved. 
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V. COMPARISON OF MODELS 


A. ANALYSIS FOR SAME DATA 

In this chapter, we continue to consider the allocation 
of spare parts to maximize the system performance in the 
different allocation models. In this thesis, we have looked 
at three models: the units short model, the time-weighted 
units short model, and the availability model. Since each 
model attempts to reduce stockouts as much as possible the 
allocations generated by the models are strongly correlated. 
This is especially true for the availability model and the 
MSRT model since availability is a function of MSRT. However, 
we saw earlier that the allocations from the models are not 
necessarily the same. 

Assume that a weapon system consists of 10 items, the 
demands are Poisson distributed and MTTR,, Cis Tis E. are 
known constants and a budget constraint of the weapon system 
is $1170. The optimal allocations for the three models are 
shown in Table IX. As can be seen, the TWUS model is more 
sensitive to the lead times than are the other two models 
(see items 5, 6, and 7). 

The units short model is more sensitive to the price of 
the item than are the other two models. For item 9 the units 
short model bought nothing, but the TWUS model and the availa- 


bility model allocated 2 and 3 items respectively. All 
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TABLE IX 


THE ALLOCATIONS OF SPARE PARTS FOR THE 
THREE DIFFERENT MODELS 


Units 
Item dj cost Ess. T; MT TR Short TWUS Avail. 
(yr) ($) x) (ea Model Model Model 
1 SO 10 il it 0.0137 3 3 4 
2 Oeed 10 ili i OOS Bien i: il 2 
3 Se 0 3 il i 0 0137 24 20 22 
4 15.0 3 3 L 0.0274 26 22 24 
5 S20 10 I 0.5 0.0274 4 5 4 
6 3as0) 5 3 0.5 0.0274 6 4 6 
7 10.0 50 il 0.2 0.0054 0 0 3 
8 Oe 0 50 1 il 0.0411 8 i 2 
g 240 50 ili ik 0301.37 0 2 3 
10 20 100 4 2 OPES 9.19) 5 > 5 


three models are highly affected by the essentiality code. 
This is illustrated by a comparison of items 9 and 10. Table 
X presents the corresponding values of the three objective 


functions. 
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TABLE X 


THE COMPARISON OF OBJECTIVE VALUES FOR THREE MODELS 


OBJ 

FN UNITS 
MODEL SHORT TWUS AVAILABILITY 
UNITS SHORT O,1523 0.0458 0.0483 
TWUS 0.1527 0.0353 0.0687 
AVAILABILITY 0.245969 O.2079)5 0.0709 


The above table was established by computing each objective 
function for the allocations determined by the three differ- 
ent procedures. Comparing the results of the three models, 
the TWUS model seems to do the best job considering all three 
objective functions. However, no general conclusions can be 
drawn about the preference of the TWUS model for other 
eetee ons. 

One needs to determine which objective function most 
closely matches a servicers' feeling about how operational 
readiness is affected by stockouts and delays in satisfying 


stockouts. 


B. DISCUSSION OF SIMILARITIES 

In the budget allocation problem there are many factors 
which affect the allocation such as demand, lead time, cost, 
time to repair and essentiality. The three models share simi- 


lar properties. First of all, as can be seen in the above 
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example, all models tend to stock the cheap, high demand 
items in favor of expensive low demands items. This is 
because of the models attempts to get the biggest benefit 
per dollar spent. Potential benefit per additional unit 
increases with an items demand rate. Second, items having 
high essentiality code are given preference, as is the in- 
tent of essentiality assignment schemes. Essentiality 
weighting is one way to counter the preference given the high 
demand low cost item observed earlier. It is frequently the 
case that the most critical items are low demand expensive 
items. Without the essentiality weighting such items would 


be neglected by the type of models examined in this thesis. 
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VI. USE OF THE MODELS FOR BUDGET DETERMINATION 


A. EFFECTIVENESS VS. BUDGET 

The models that we have discussed have attempted to opti- 
mize performance subject to a budget constraint. We have 
assumed that the budget was given. There are many ways in 
which budgets are determined. However, budgeting people and 
inventory managers alike often express the desire to have a 
methodology that they can use to determine the amount of 
money that should be provided. 

In most cases the amount is determined historically by 
giving an amount equal to what has been provided in the past 
for similar systems or perhaps by giving a little more or 
less based on judgement or financial constraints. There is, 
however, a strong interest brought about by Congressional 
pressures to relate resources to readiness. Congress wants 
to know "how much money is needed to support our weapon sys- 
tems at a specified level of performance." In this chapter we 
show how the models developed earlier in this thesis can be 
used in just this manner. 

Specifically, we show how the models that we have developed 
can be modified easily to determine the minimum budget re- 
ase to provide a specified level of logistics performance. 

The models developed earlier can each be run for a range 


of budget levels producing for each given budget an allocation 
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and a predicted overall level of performance. Figure 3 illus- 
trates this for the case in which the performance measure is 
pseudo~-availability. As expected, the curve shows that 
availability is a non-decreasing function of budget with 
decreasing marginal returns. This can be done also for the 
time-weighted units short model or any of the other models 
discussed in this thesis. In all cases we would obtain a 
similar display. Performance is a monotonic function of 
budget with decreasing marginal returns. 

Figure 4 displays a similar result for the case in which 
the performance measure is MSRT. Each point on the curve 
represents an optimal level of performance for a given budget. 
For this example displayed, Figure 3, there is a little bene- 
fit to be gained by increasing the budget above $2500. How- 
ever there is a dramatic increase in effectiveness obtained 
by increasing the budget from $1000 to $2000. This 1s pre- 
cisely the sort of information needed to make intelligent 
budgeting decisions. Of course some decision maker must 
decide if the increase in effectiveness is worth the addi- 
tional expenditure. 

If a specific level of effectiveness is specified, one 
can graphically determine the amount of budget required by 
simply moving horizontally across the graph from the speci- 
fied level of effectiveness until the curve is intersected 


and then down to the budget axis. 
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AVAILABILITY VS BUDGET 


6.350 


040 


as 


030 


AVAIL. 
& 


0.20 
0.8 
0.99- 
6.05 


0.00 = 
© S00 100 S00 2000 71300 3000 3500 4000 4500 3000 
BUDGET 


Figure 3. Total Availability Vs. Budget Curve 
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The next section determines analytically the minimum 
amount budget required by solving a companion problem to 


the problems discussed earlier in this thesis. 


B. COMPANION PROBLEMS 

In the previous chapters we have concentrated on the 
optimization of system effectiveness with a budget constraint. 
For many weapon systems such as air detection radars, missiles 
and nuclear delivery systems, the system performance is so 
important that the necessary budget will be provided to 
attain whatever performance is deemed necessary. 

For such systems it is reasonable to restate the optimi- 
zation problem to determine the minimum budget required to 
satisfy a specified level of performance. Consider, for 
example, the availability optimization problem and the com- 


panion problem: 


n 
(D1) min a ers 


n 
S.t. I A, > fee es GN 


where L is the minimum performance level for a weapon system. 


For the MSRT model case the corresponding problem is: 


(D2) min ) c.s, 


n 
; PBL ex eat 
ices ). MSRT,(s,) < R i= 1 n 


i=l 


a 





where R is the maximum allowable cumulative supply response 
time for the weapon system. 

Problems (Dl) and (D2) can be solved using the same 
methods explained in Chapters III and IV. In the above 
models the budget is determined so that the system require- 
ment for availability or main supply response time can be 
achieved. 


For problem (D2) the total cost 1S minimized when 


n 
)ORMSRT(s.) = R. 
ee oe 
1=1 
Sometimes a minimum allowable supply response time is 
required for each item. In such a case multiple constraints 


could be specified. This is illustrated below: 


n 
(D3) min af, c.s, 
MSRT, (sj) < Ry 
MSRT. (s,) < Ro 
MSRT_(s_) < RE 


where R, is the maximum allowable supply response time, 


Ie oy ss JN 
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To solve (D3), find the smallest S; such that 


MSRT, (s,) = Ry 
MSRT.(s.) = Ro 
MSRT_(s_) = R, 


This problem is solved easily using the same procedures 
which we discussed. 

So far we have discussed many different ways to apply 
the theoretical models to practical use of models for budget 
determination decisions. There is no unique method which 
gives us an optimal result. So the user of these models 
should choose one of possible methods so as to maximize the 


system performance or minimize the total cost. 
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VII. CONCLUSIONS 


It is concluded that the various measures of effectiveness 
can be used in the budget constrained multi-item inventory 
system with stochastic demands. We have examined some of 
the more reasonable measures like minimization of units short, 
minimization of time-weighted units short and maximization 
of system availability. We have also looked at models which 
incorporate essentiality weights into each of the models. 

In order to solve budget allocation stockage problems a 
feasible, efficient method of effecting line item inventory 
control is available using an adaptation of Everett's 
Generalized Lagrangian Multiplier method. Further, the use of 
a G.L.M. procedure provides valuable information for system 
Managers as to relative effectiveness of additional procure- 
ment funds, versus additional transition processing capability. 
The final value of Lagrangian multipliers can be interpreted 
as the amount of improvement of the objective function per 
unit dollar spent. 

The models discussed in this thesis are all more likely 
to stock cheap, high demand items than expensive, low demand 
items. Such is the nature of budget constrained optimization 
problems. If a system manager wishes to maintain enough 
stock for an item having low demand, high cost, his only 


alternative in our models is to assign a high essentiality 
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code for the item. The essentiality code has the effect of 
reducing the ratio C/E as opposed to C. In'the solution 
procedure for each model, the assigned essentiality code 
directly affects the allocation for the item. 

_We have shown how the models can be used as a tool to 
determine the amount of budget. A simple graphical procedure 
allows a decision maker to determine the minimum budget re- 
quired to search a specified level of performance. The opti- 
mization model is run several times to generate a plot of 
performance vs. budget. Each point of the curve represents 
the effectiveness for the optimal allocation of a given budget. 
A manager can, first, determine an appropriate system per- 
formance level and read from the curve the budget required 
to achieve the effectiveness. 

Further analysis to improve these models may be possible. 
For instance, it would be useful to have an automatic search 
algorithm for the Lagrangian multipliers for a multiple 
constrained problem. It may also be possible to relax the 
assumptions for a constant lead time or a constant mean time 
to repair. These single period inventory may expand to time- 
dependent multi-item, multi-echelon, multi-indenture inventory 


systems. 
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APPENDIX A 


COMPUTER PROGRAM FOR INTERACTIVE SEARCH 
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APPENDIX B 


COMPUTER PROGRAM FOR AUTOMATING SEARCH 
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COMPUTER PROGRAM FOR MARGINAL ANALYSIS 
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